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Abstract 

A recent analysis of the 4 He abundance determined from observations of extragalactic HII 
regions indicates a significantly greater uncertainty for the 4 He mass fraction. The derived 
value is now in line with predictions from big bang nucleosynthesis when the baryon density 
determined by WMAP is assumed. Based on this new analysis of 4 He, we derive constraints 
on a host of particle properties which include: limits on the number of relativistic species 
at the time of BBN (commonly taken to be the limit on neutrino flavors), limits on the 
variations of fundamental couplings such as a em and Gn, and limits on decaying particles. 



1 Introduction 



Big bang nucleosynthesis (BBN) is one of the most sensitive available probes of physics 
beyond the standard model. The concordance between the observation-based determinations 
of the light element abundances of D, 3 He, 4 He, and 7 Li [1], and their theoretically predicted 
abundances reflects the overall success of the standard big bang cosmology. Many departures 
from the standard model are likely to upset this agreement, and are tightly constrained [2]. 

The 4 He abundance, in particular, has often been used as a sensitive probe of new physics. 
This is due to the fact that nearly all available neutrons at the time of BBN end up in 4 He and 
the neutron-to-proton ratio is very sensitive to the competition between the weak interaction 
rate and the expansion rate. For example, a bound on the number of relativistic degrees 
of freedom (at the time of BBN), commonly known as the limit on neutrino flavors, N u , 
is derived through its effect on the expansion rate, H oc ^fgl [3]. However, because the 
calculated 4 He abundance increases monotonically with baryon density (parameterized by 
the baryon-to-photon ratio, r\ = rib/n^), a meaningful limit on iVj, requires both a lower 
bound to i] and an upper bound to the primordial 4 He mass fraction, Y p [4]. Indeed, for a 
fixed upper limit to Y p , the upper limit to N v can be a sensitive function of the lower limit 
to i], particularly if rj is small [4,5]. 

The recent all-sky, high-precision measurement of microwave background anisotropies 
by WMAP [6] has opened the possibility for new precision analyses of BBN. Among the 
cosmological parameters determined by WMAP, the baryon density has been derived with 
unprecedented precision. The WMAP best fit assuming a varying spectral index is f^/i 2 = 
0.0224 ± 0.0009 which is equivalent to ^io.cmb = 6.14 ± 0.25, where 7710 = 10 1( V This 
result is sensitive mostly to WMAP alone but does include CMB data on smaller angular 
scales [7] , Lyman a forest data, and 2dF redshift survey data [8] on large angular scales. This 
result is very similar to the corresponding value obtained from combining WMAP with SDSS 
data and other CMB measurements, which gives VL^h 2 = 0.0228io;ooJg [9] and corresponds 
to r]i = 6.251q!22- Using the WMAP data to fix the baryon density, one can make quite 
accurate predictions for the light element abundances [10-13]. At the WMAP value for 77, 
the 4 He abundance is predicted to be [13]: 

Y p = 0.2485 ± 0.0005 (1) 

On the other hand, accurate 4 He abundances have been and continue to be difficult to 
obtain. It is recognized that there are many potential sources of systematic uncertainties 
in the derived 4 He abundance [14, 15]. As a result there exists a wide range of derived 
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primordial 4 He abundances which have typically been relatively low compared with ((TJ. 
Recently, a reanalysis [16] of the 4 He data [17,18] has led to a significant enlargement in the 
statistical uncertainty as well as a potential shift in the mean value. A representative result 
of that analysis is 

Y p = 0.2495 ± 0.0092 (2) 

Conservatively, it would be difficult to exclude any value of Y p inside the range 0.232 - 0.258. 

Because much of the previous work was based on relatively low values of Y p , tension 
between the value of 77 inferred from either D/H or WMAP, and 4 He gave rise to very tight 
constraints on N v and on other particle properties. In light of the newly suggested range for 
Y p [16], it is important to reexamine the constraints on physics beyond the standard model. 
Potential limits from D/H have been discussed recently [10,12,13], and we will just quote 
those results below in comparison with the results derived here. At present, it is not possible 
to use 7 Li to obtain constraints. This is due to 1) the large uncertainty in the BBN prediction 
of the 7 Li abundance, and 2) to the current discrepancy between the BBN prediction and 
the observational determination of the 7 Li abundance (see e.g. [10-13,19]). 

2 The 4 He Abundance 

The 4 He abundance has had a somewhat checkered history over the last decade. Of the 
modern determinations, the work of Pagel et al. [20] established the analysis techniques that 
others were soon to follow [21]. Their value of Y p = 0.228 ± 0.005 was significantly lower 
than that of a sample of 45 low metallicity HII regions, observed and analyzed in a uniform 
manner [17], with a derived value of Y p = 0.244 ± 0.002. An analysis based on the combined 
available data as well as unpublished data yielded an intermediate value of 0.238 ± 0.002 
with an estimated systematic uncertainty of 0.005 [22]. An extended data set including 89 
HII regions obtained Y v = 0.2429 ± 0.0009 [18]. However, the recommended value is based 
on the much smaller subset of 7 HII regions, finding Y p = 0.2421 ± 0.0021. As seen in table 
6 of [18], changing the assumed value of the equivalent width of He absorption for one of 
the observed He wavelengths by 0.1 A changes the derived abundance significantly to Y p = 
0.2444 ± 0.0020. This change of AY P = 0.0023, is indicative of the importance of systematic 
errors. 

4 He abundance determinations depend on a number of physical parameters associated 
with the HII region in addition to the overall intensity of the He emission line. These include 
the temperature, electron density, optical depth and degree of underlying absorption. A self- 
consistent analysis may use multiple 4 He emission lines to determine the He abundance, 
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the electron density and the optical depth. In [17], five He lines were used, underlying He 
absorption was assumed to be negligible and temperatures based on OIII observations were 
used. 

A very accurate helium abundance for the HII region NGC 346 in the Small Magellanic 
Cloud was derived with a value of Y p = 0.2345 ± 0.0026 [23]. Knowing that the OIII 
temperatures are systematically high, they use the He I emission lines to solve for the electron 
temperature. Recently, the spectra of five metal poor HII regions - NGC 346 and four 
regions reported in [17] have been reanalyzed [24]. After considering the effects of additional 
physical processes (e.g., collisional excitation of the Balmer lines), a higher determination of 
Y p = 0.239 ± 0.002 was found. 

The question of systematic uncertainties was addressed in some detail in [14]. It was 
shown that there exist severe degeneracies inherent in the self-consistent method, particularly 
when the effects of underlying absorption are taken into account. A sixth He line was 
proposed to test for the presence of underlying He absorption. However, even in the six-line 
method, one can not escape the degeneracies present in the solutions. In particular, solutions 
with no absorption and high density are often indistinguishable (i.e., in a statistical sense 
they are equally well represented by the data) from solutions with underlying absorption 
and a lower density. In the latter case, the He abundance is systematically higher. These 
degeneracies are markedly apparent when the data is analyzed using Monte-Carlo methods 
which generate statistically viable representations of the observations. When this is done, 
not only are the He abundances found to be higher, but the uncertainties are also found to 
be significantly larger than in a direct self-consistent approach. 

In [16], the Monte-Carlo method was applied to seven of the highest quality observations 
from the sample in [17]. As expected, systematically higher He abundances were found, 
with significantly larger uncertainties. The results of a regression (to zero metallicity) led 
to the primordial He abundance given in Eq. (j2J). With this value of Y p (particularly with 
the stated uncertainty), there is clearly no discrepancy between He observations and BBN 
predictions of Y p at the WMAP value for rj. It was stressed however, that although this 
method found higher He abundances, one could not exclude the lower abundances found by 
other methods. 

3 Standard BBN 

Key to BBN analysis is an accurate determination of BBN theory uncertainties, which are 
dominated by the errors in nuclear cross section data. To this end, several groups have 
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determined reaction rate representations and uncertainties. Smith, Kawano and Malaney [25] 
presented the first detailed error budget for BBN, generally assuming constant relative errors. 
In more recent work [26] , uncertainties were propagated based on available nuclear data into 
the light element predictions. The NACRE collaboration presented a larger focus nuclear 
compilation [27], meant to update the previous astrophysical standard [28]. However, their 
"high" and "low" limits are not defined rigorously as 1 or 2 sigma limits (see [29,30] for its 
impact on BBN). In an attempt to increase the rigor of the NACRE errors, we reanalyzed 
[29] the data using NACRE cross section fits defining a "sample variance" which takes into 
account systematic differences between data sets. 

Since then, new data and techniques have become available, motivating new compila- 
tions. Within the last year, several new BBN compilations have been presented [11-13], 
the latter being one of the most rigorous and exhaustive efforts to determine reliable rate 
representations and meaningful uncertainties and we adopt this compilation here. 

To compare theory with observations, we will adopt the 4 He results discussed in £0 For 
brevity, we will simply adopt the D observations discussed in [31]. We use the D abundance 
constraints based on the best five measurements of D/H in QSO absorption line systems 

D/H A = (2.78 ± 0.29) x 10" 5 (3) 

BBN is tested with the comparison between light element abundance observations, BBN 
theory predictions, and the subsequent allowed ranges of baryon density, which has now been 
independently measured by CMB anisotropy experiments [10]. One way to perform this test 
is to use the CMB range of r\ as an input to the BBN calculation, and to compute the resulting 
ranges in the light elements. This procedure is illustrated in Figure ^ As is well known, the 
agreement with D is excellent. The 4 He overlap is also essentially perfect, though the larger 
uncertainty in 4 He renders this agreement a less powerful test, and the near-coincidence of 
central values is fortuitous. Nevertheless, the new 4 He analyses does bring this nuclide in 
good agreement with the CMB and D. This consistency was not guaranteed and marks a 
success of BBN and of cosmology. 

4 Beyond the Standard Model 

For several cases of interest, it will be useful to define a dimensionless cosmic "speed-up" 
factor £ = if new /if st d, where H = a/a is the Hubble expansion rate; £ = 1 then represents 
the unperturbed case. The expansion rate itself is given by the Friedmann equation, which 
for a flat universe is H 2 = (87r/3) Gj^p, where p is the total mass-energy density. Thus the 
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Figure 1: Likelihood distributions for light element abundances. The dark (blue) shaded 
regions are the BBN predictions given the CMB-determined r\ values. The light (yellow) 
shaded regions correspond to the light element observations: the new Y p and world average 
of D/H^ values of eq. © are represented by the yellow regions. 



speed-up factor evolves as £ = y (G np) new /(GNP)std- For the case of a radiation-dominated 
universe, we have p oc g*T 4 , where g>* = 2 + 7/2 + 7N u /4 counts the relativistic degrees of 
freedom in photons, e ± pairs, and N v neutrino species. 

4.1 Constraints on N v 

We first consider the canonical extension of standard BBN, in which there are N v effectively 
massless (m u «C 1 MeV) left-handed neutrino species. The increase in the speed-up factor is 



where 6N V = N u — 3. This in turn changes the weak freezeout temperature and ultimately 
affects all of the light elements. Until recently, the effect on the 4 He abundance was the only 
measurable consequence, but D/H measurements are now sufficiently accurate that D/H also 
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Figure 2: The likelihood distribution for N v based on the WMAP value of 7710 = 6.14 ±0.25 
and Y p from Eq. H (shaded), WMAP and D/H (dashed), WMAP and both Y p and D/H A 
(dotted). We also show the result without the imposing the WMAP value for 77 (long dashed). 

has an important sensitivity to N u [32]. 

For a fixed value of 7710 = 6.14 ± 0.25 and the He abundance given in (J2J), we show the 
likelihood distribution for N v by the shaded region in Fig. |21 Also shown for comparison are 
the likelihood distribution based the WMAP value of 77 using D/H alone, Y p and D/H, and 
the result based on BBN alone. Despite the increased uncertainty in the He abundance, it 
still provides the strongest constraint on N u . D/H is nonetheless becoming competitive in 
its ability to set limits on N v . 

Figure shows the joint limits on 77 and N v based on D and 4 He. We see that the 
4 He contours are nearly horizontal, which arises from the weak (logarithmic) sensitivity of 
Y p to 77, as opposed to a stronger, linear sensitivity to 5N U . Thus 4 He by itself is a poor 
baryometer but an excellent probe of nonstandard physics. On the other hand, the D/H 
contours have a steep slope, indicating a strong sensitivity to 77 which is the origin of the 
D/H power as a baryometer. The non- vertical nature of the slope does however indicate a 
correlation between the D/H sensitivity to 77 and N v . Thus by combining D and 4 He we 
can expect to arrive at strong constraints on both parameters. Numerical results appear in 
Tabled where we see that these light elements alone constrain 77 to within about 10%, and 
fix N u to within about 20%, both at la. Note that the contour ellipses in Figure El have a 
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Figure 3: a) BBN-only constraints on rj and N u . The thickest (thinnest) curves correspond 
to la (3a) limits. The nearly vertical (blue) curves are limits due to D/H, nearly horizontal 
(red) curves are for 4 He, and the closed (black) contours combine both, b) As in a), with 
the CMB 7] information included. 

slight positive tilt, corresponding to a small positive correlation between rj and N u . 

Also appearing in Tabled we have shown the 95% upper confidence limits placed on the 
effective neutrino number, 5N UtTnax , assuming that N v > 3.0 or 5N U > 0.0. The constraints 
presented suggest a robust upper bound of 1.6 with 95% confidence. We next introduce the 
CMB information on rj] this tests the overall consistency, but as we have already shown, the 
agreement is good for the standard N v = 3 case. Note that CMB anisotropics also have some 
sensitivity to N u , though this is at the moment significantly weaker than the light element 
sensitivity. We do not use this additional information, which would slightly strengthen the 
constraints on N v , but would not affect the rj limits (where the CMB impact is largest) due 
to the independence of the CMB limits on 77 and N v [33]. 

Figure Eb shows the impact of the CMB on the rj and N v constraints. We see that 
the dominant effect is that the CMB narrows and steepens the combined contours; this 
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Table 1: The table shows constraints placed on N v and r\ by various combinations of observa- 
tions. Shown are the 68% confidence limits determined by marginalizing the 2-D likelihood 
distribution C(t),N u ). Also shown are the 95% upper limits on 8N V = N u — 3, given that 
5N V > 0. 
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reflects the very tight CMB constraint on 77. Table ^ shows the impact of the CMB on the 
rj and N v constraints. The resulting precision on r\ is roughly doubled, to about a 4% (!) 
measurement, dominated by the CMB contribution but for which the D/H contribution is 
not negligible. The precision of the N u constraint remains essentially the same, reflecting 
both the dominance of Y p in determining N u , as well as the near-independence of Y p on rj. 

In [16], it was noted that the primordial value of the 4 He abundance based on a regression 
with respect to O /H was only marginally statistically more significant that a weighted mean 
which yields Y p = 0.252 ± 0.003. This result is also obtained using a Bayesian analysis in 
which the sole prior is the increase in 4 He in time [34] . The combination of the 4 He abundance 
based on the mean value and the CMB value for rj gives N u = 3.27 ± 0.24 with a 95% upper 
limit 5N ViTnax = 0.7. Recall, the constraint on SN U)Tnax assumes N u > 3 or ^A"^ > 0. 

In all cases the preferred values for A^ are consistent with N u = 3, and in many cases 
are much closer to N v = 3 than la. This restates the overall consistency among standard 
BBN theory, D and 4 He observations, and CMB anisotropics. It also constrains departures 
from this scenario. Our combined limit using BBN ± light elements + CMB limit is: 

2.67 <N U < 3.85 (5) 

at 68% CL. 

4.2 Constraints on the Variation of Fundamental Constants 

As noted earlier, BBN also placed interesting limits on possible variations of fundamental 
constants. Indeed, almost every fundamental parameter can be constrained by BBN if 
it affects either the expansion rate of the Universe, the weak interaction rates prior to 
nucleosynthesis, or of course the nuclear rates themselves. As quantitative examples of the 



8 



constraints which can be derived, we focus here solely on variations of Newton's constant 
and the fine-structure constant. Here, we simply note that many other constraints have been 
considered in the past (for a recent review see: [35]). 

4.2.1 Newton's Constant, G N 

Strictly speaking, it makes no sense to consider the variation of a dimensionful constant such 
as G N (see e.g. [36]). We include the discussion here for the purposes of comparison with 
previous constraints. We also note that the limit set here could be interpreted as a limit on 
the gravitational coupling between two protons {G^m^) in a framework where we have chosen 
m p (and all other particle masses) to be constant. Early constraints [37] on Gn relied mainly 
on 4 He abundance observations. Recently, the D/H abundance was used [38] in conjunction 
with the WMAP determination of rj to set a limit on AGn/Gn,o of about 20% from the 
time of BBN. Assuming a simple power law dependence Gn ~ t~ x , x was constrained to the 
range -0.004 < x < 0.005 implying -4 x 10~ 13 yr _1 < G Nfi /G Nfi < 3 x lO^yr -1 The use 
of D/H was motivated in part by the previously discrepant value for Y p . 

Although the uncertainty in the 4 He abundance has been argued to be significantly larger 
than past values (0.009 vs 0.002) [16], the resulting bounds on G N are still interesting. The 
limit on 5N V of —0.60 < 5N V < 0.82 from the Y p + CMB constraint, can be translated 
directly to a bound on the speed-up factor: 0.949 < £ < 1.062. Any variation in Gn can be 
expressed through £ 

- 0.10 <^=£ 2 -l< 0.13 (6) 

or perhaps more simply put, the limit on the variation in Gn can always be related directly 
to the limit on N v through 

Gn,o 43 

If one makes the common assumption that Gn ~ t~ x , we obtain, —0.0029 < x < 0.0032 
and -2.4 x lO^yr" 1 < G Nfi /G Nfi < 2.1 x lO^yr" 1 (using t = 13.7 ± 0.2 Gyr [6], and 
tBBN ~ 100 sec). Thus despite the increased uncertainty in Y p , the 4 He abundance still 
provides the strongest possible constraint on Gn- 

4.2.2 The Fine-Structure Constant, a 

There has been a great deal of activity surrounding possible variations of the fine-structure 
constant, motivated largely by a reported observational analysis of quasar absorption systems 
which has been interpreted as a variation in a [39]. We note that other observations using 
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similar methods [40] have not confirmed the variation in a, and other interpretations based 
on the nucleosynthesis of heavy Mg isotopes in the absorbers may also explain the data [41]. 
Other constraints from the CMB [42], the Oklo reactor [43,44], and meteoritic abundances 
[44,45] have also been derived. Once again, our goal here is to update the BBN bound on 
variations in a using the newly derived value of Y p . 

If we assume that only a is allowed to vary (i.e., we assume that all other fundamental 
parameters are held fixed), the dominant contribution to a change in Y p comes from the 
variation in the neutron-proton mass difference, Q = m n — m p [46]. The 4 He abundance 
can be estimated simply from the ratio of the neutron-to-proton number densities, n/p, by 
assuming that essentially all free neutrons are incorporated into 4 He. The neutron-to-proton 
ratio at weak freezeout is (n/p)f ~ e~Q' T f , modulo free neutron decay, where Tf is the 
temperature at which the weak interaction rate for interconverting neutrons and protons 
falls below the expansion rate of the Universe. 

Variations in Q leads to a variation in Y p given approximately by 

AF ~ AQ (R) 
— ~ ~Q (8) 

One can write the nucleon mass difference as 

Q ~ aaA QC D + bv (9) 

where a and b are dimensionless constants giving the relative contributions from the electro- 
magnetic and weak interactions. In (|9*)l. v is the standard model Higgs expectation value. A 
discussion on the contributions to Q can be found in [47] . The constants a and b are chosen 
so that at present the two terms contribute -0.8 MeV and 2.1 MeV respectively. Eqs. (|8ji 
and © can be combined to give 

^ K 0.6^ (10) 
Y a 

Thus the current uncertainty in the observational determined value of Y p leads to a bound of 
\Aa/a\ < 0.06. If changes in a are correlated to changes in other gauge or Yukawa couplings, 
this limit improves (in a model dependent way) by about 2 orders of magnitude [48] . 



4.3 Limits on Decaying Particles 

An exotic scenario often considered is that of late-decaying particles (j~x ~ 10 8 sec) [49]. 
The particles are assumed to decay electromagnetically, meaning that the decays inject elec- 
tromagnetic radiation into the early universe. If the decaying particle is abundant enough 
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or massive enough, the injection of electromagnetic radiation can photo-erode the light ele- 
ments created during primordial nucleosynthesis. The theories we have in mind are generally 
super symmetric, in which the gravitino and neutralino are the next-to-lightest and lightest 
super symmetric particles, respectively, but the constraints hold for any decay producing 
electromagnetic radiation. We thus constrain the abundance of such a particle given its 
mean lifetime tx- The abundance is constrained through the parameter (x = 2E in jnx/n^. 
We can see there is a degeneracy between the relative abundance nx/n 7 and the injected 
energy E in j. Specific theories can relate the lifetime with the mass of the decaying particle 
and thus the injected energy E in j, however, we will restrict ourselves to the general case. 

The constraint placed by the 4 He abundance comes from its lower limit, as this scenario 
destroys 4 He. The maximum value our abundance parameter can take is proportional to 
the 4 He constraint, Cx,max( 4 He) oc (Y BBN - Y min )/Y BBN . Here Y BBN is the predicted 4 He 
abundance given r\ and Y min is the minimum allowed value from observations. Using Y min = 
0.232, we find: 

C*( 4 He) < 2.1 x 10- 10 GeV (|^) (j^)^ for r x > 10* sec (11) 

However, in this scenario the limits on the tertiary production of 6 Li provide stronger con- 
straints on late-decaying particles, with (x < 5.1 x 10~ 12 GeV following the same scalings 
with 7710 and Tx as 4 He. 

4.4 Other Bounds 

Models of neutrino masses with new right-handed interactions are also subject to the con- 
straint from 4 He. If we assume right-handed neutrinos are present and are light enough 
(< 1 MeV) to count as extra relativistic degrees of freedom during BBN, they are subject to 
the constraint on 5N U < 1.60. Assuming there are 3 standard model neutrinos and 3 right- 
handed neutrinos, we can relate the limit on 5N U to the temperature ratio of right-handed 
to left-handed neutrino temperatures: 6N V = ?>{T VR /T VL ) i [50]. We find T Vr /T Vl < 0.85 with 
95% confidence. In order to change this ratio from unity, right-handed neutrinos must decou- 
ple before epochs of entropy release, generally due to particle annihilation. To accommodate 
the constraint on the temperature ratio, z/r must decouple at least before the annihilation 
epoch of pions, if not before the quark-hadron transition, with T dec > 140 MeV. 

As a consequence, right-handed interactions governed by some new mass scale M R , can 
also be constrained. Requiring that these interactions decouple before the pion annihilation 
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epoch and recalling that Td ec scales as M^, we find that: 

/ Ta \ 3/4 

M »'« >3 - 3Tev (i4o!y <i2 > 

Note this constraint is stronger than limits currently set by high-precision electroweak con- 
straints (M Wr > 715 GeV [51]). For recent analyses of models of this type, see [52]. 

BBN can also constrain the presence of vaccum energy during nucleosynthesis, via the 
expansion speedup provided by the additional energy density. The presence of vaccuum 
energy during BBN is motivated both by considerations of dark energy [53, 54] as well as 
inflation. The new vacuum energy component is typically modelled as some scalar field 
cf). The evolution of the scalar field is controlled by its potential as well as the Friedmann 
equation. However, for an important class of models, there is a "tracker" solution which 
is an attractor, and leads the vaccum energy density to be proportional to the dominant 
radiation density [53]. In this case, the vacuum energy acts essentially as an extra neutrino 
species, with 

0,(1 MeV) = — = 76N j A , (13) 

^ } Prad + P^ 10.75 + 75^/4 1 ; 

Using our limit 5N V < 1.60, we find that f^(l MeV) < 0.21. For the case of an exponential 

potential V(<f>) = M^ x e~ x ^ Mpi , this constraints the coupling A = 2/y^ > 4.4. 

5 Summary 

A new and detailed assessment [16] of the observed primordial 4 He abundance, and its 
uncertainties, has important implications for cosmology and BBN generally, and for early 
universe and particle physics in particular. The observed 4 He abundance is now found to 
be consistent with the r\ value given by D, leaving 7 Li alone in discordance. Moreover, both 
D and now 4 He are consistent with the precision r\ range determined by recent observations 
of CMB anisotropies. The newfound 4 He agreement arises primarily because the new and 
more detailed error budget is also larger than previous estimates. Nevertheless, with 4 He 
in concordance, it now rejoins D and the CMB as a probe of physics beyond the Standard 
Model. 

We have surveyed the impact of the new 4 He analysis on nonstandard physics. The classic 
neutrino counting argument is found to give relaxed though non-trivial limits to 5N U (eq. El 
and Table Q). We find that D and 4 He observations each make similar contributions to these 
limits. Thus, if one of these observations can be signficantly improved, it will dominate our 
ability to probe exotica. 



12 



These limits immediately translate to tighter constraints on variation of the fundamental 
constants. They also constrain the presence of vaccuum energy (in the "tracker" regime) 
during BBN. Finally, we note that the new 4 He limit on decaying particle sceanrios remains 
weaker than the constraints placed by 6 Li; with a factor ~ 2 improvement in precision, 4 He 
would be competetive. 

We thus urge continued effort to (realistically!) improve the precision of the observations 
of primordial abundances. A sharpening in either D or 4 He will immediately tighten the 
limits we place here. And of course, the outstanding problem with 7 Li continues to demand 
explanation. 
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